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1 Abstract 

This paper presents the major characteristics of 
a simulation package capable of performing a 
complete narrow and wideband analysis of the 
mobile satellite communication channel in ur- 
ban environments, for any given orbital config- 
uration. For the RF frequency range the model 
has been designed to be applicable (1 up to 60 
GHz), the wavelength-to-average urban geomet- 
rical dimension ratio has required the use of the 
Geometrical Theory of Diffraction (GTD) [1], 
extended to include effects of non-perfect con- 
ductivity and surface roughness. Taking ad- 
vantage of the inherent capabilities of such a 
high frequency method, we are able to provide a 
complete description of the electromagnetic field 
at the mobile terminal. Using the information 
made available at the ray-tracer and GTD solver 
outputs, the Land Mobile Satellite (LMS) ur- 
ban model can also give a detailed description 
of the communication channel in terms of power 
delay profiles, Doppler spectra, channel scatter- 
ing functions and so forth. Statistical data, e.g. 
cumulative distribution functions, level crossing 
rates or distributions of fades are provided too. 
The user can access the simulation tool through 
a Design-CAD user-friendly interface by means 


of which he can effectively design his own urban 
layout and run consequently all the envisaged 
routines. The software is optimised in its execu- 
tion time so that numerous runs can be achieved 
in a considerable short time. 

2 Introduction 

Urban areas will likely represent a significant 
market for the conventional mobile communi- 
cation services and for the new concepts of per- 
sonal communications through hand-held termi- 
nals. The propagation of an electromagnetic 
field in such environments is, on the other hand, 
a very complex phenomenon to simulate. To 
date the analyses of the link impairments and 
the estimation of their impact on the LMS sys- 
tem performance have been mainly based upon 
experimental data and empirical models, [2]-[4]. 
The development of a deterministic model not 
directly related to specific environmental urban 
scenarios and designed on canonical electromag- 
netic laws is therefore strongly needed. 

Based on this sort of considerations and hav- 
ing in mind to develop a user-friendly predic- 
tion tool to be used not only by propagation 
engineers but also by LMS system planners, the 
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European Space Agency (ESA) decided to place 
a study contract with two Italian companies, 
Space Engineering and Ingegneria dei Sistemi 
(IDS), in 1992, [5], The main features of the 
LMS prediction tool developed under this con- 
tract are hereafter described. 

3 The LMS urban predic- 
tion tool 

The GTD-based prediction tool for LMS sys- 
tems in urban environments can be considered 
as the result of the interaction of several func- 
tional blocks, each of them properly designed in 
order to save computational time and increase 
modularity and transportability. The core of 
the simulation package is represented by the Ray 
Tracer and the GTD Solver codes. In addition 
we also have an Urban Area Modeler, an Elec- 
tromagnetic (E.M.) Mesher and a Post Proces- 
sor Unit. A functional block diagram of the 
package is reported in Fig. 1 while a detailed 
description at subsystem level is given in the 
following sections. 

3.1 The Ray Tracer 

In most of the existing GTD prediction pack- 
ages (e.g. NEC-BSC) ray tracing and electro- 
magnetic field computation are performed in 
the same logical step; the ray tracing parame- 
ters calculated at any observation point are not 
saved for further processing (e.g. at a differ- 
ent RF frequency) and this inevitably results in 
a great computational inefficiency particularly 
in complex environments such as urban areas. 
In our LMS prediction tool, the two operations 
are kept separate and are performed by two dis- 
tinct modules. All the ray tracing parameters 
are saved in a file before actually starting the 
e.m. calculation. This software architecture im- 
plies a sensible reduction of the CPU time and is 
naturally suitable for optimization and upgrad- 
ability. 


The input parameters for the Ray Tracer are: 
the satellite position, the urban area model, the 
vehicle speed and path and the number of elec- 
tromagnetic interactions to be considered. The 
model has in-built the opportunity to gener- 
ate orbital configurations other than the con- 
ventional geostationary one; this feature is ex- 
tremely important to simulate a large variety of 
LMS constellations, from Low to Medium and 
Helliptical Earth Orbiting multisatellite systems 
(LEO, MEO, HEO). The urban geometry of the 
vehicle path is given through a set of straight 
lines with a user- assigned speed for each of them. 
As for the number of the electromagnetic inter- 
actions, the LMS urban prediction tool has been 
to date designed to include up to the second or- 
der, i.e. double reflected or diffracted and all 
the different combinations with the direct and 
single scattered rays. 



Fig. 1 LMS prediction tool block diagram 
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The software architecture has been anyhow 
conceived to be easily upgradable for any higher 
order of contributions. 

In order to reduce the computational time, 
the Ray Tracer considers as active part of the 
user-defined urban layout the boolean combina- 
tion of the portions contemporarily seen by the 
satellite and the mobile terminal. Once all the 
ray paths relevant to a certain observation point 
are computed, they are stored on a high capac- 
ity hard disk and the mobile position is updated. 
In Fig. 2 a sketch of the Ray Tracer active area 
for a typical urban area is reported where all 
the relevant rays (reflected, diffracted, etc.) are 
clearly visible. 



Fig. 2 Typical Ray Tracer output 

3.2 The GTD solver 

This unit performs the actual electromagnetic 
computations according to the uniform theory 
of diffraction and using, where applicable, the 
geometrical optics and the Fresnel coefficients 
for the reflection ray paths. The extension to 
non perfectly conductive materials and surface 
roughness is achieved through heuristic formu- 
lations and user-defined look-up tables for the 
complex reflection coefficients. The effect of 
leafage absorption is also modeled directly by 


the user introducing weight functions or empir- 
ical relations. 

The input parameters for the GTD Solver 
are, for any given satellite position and obser- 
vation point: the ray path history file, the e.m. 
urban model and auxiliary information, e.g. RF 
frequency. The direct, reflected, diffracted and 
refracted field components are then computed in 
amplitude and phase. This information is stored 
in a file and handed over to the Post Processor 
Unit. 

3.3 The Urban Area Modeler 
and E.M. Mesher 

The urban layout is characterised by a set of 
large and small objects representing urban pa- 
rameters, such as buildings, tunnels, overpasses, 
trees, street lamps, phone booths, parked cars, 
etc. Through the Urban Area Modeler imple- 
mented in Design-CAD, the user is able to cre- 
ate a fully controlled urban environment placing 
building by building; it is envisaged, for future 
developments, also an automatic, generation pro- 
cedure based on a set of statistical parameters 
(urbanization factors). The urban layout pre- 
viously considered for the computation of the 
active Ray Tracer area is reported in Fig. 3. 



Fig. 3 Urban Modeler output 
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The Urban Modeler has also the possibility of 
accepting as input file the layouts of actual cities 
(Paris, London, Rome, etc.) if given in forms of 
digital terrain data bases. 

Large objects are defined mainly by boxes 
and cylinders of predefined dimensions; more 
complex urban objects can be built using flat 
panels. Trees are characterised by a spheri- 
cal leafage and a cylindrical trunk while other 
small objects (nearby vehicles, dust bins, phone 
booths, etc.) are accounted for with spheres of 
given radius, scattering isotropically. All the ob- 
jects materials are defined with their complex 
permeability and permittivity. 

The E.M. Mesher is the unit interfacing the 
urban layout modeler with the electromagnetic 
solver. It performs all the required verifications 
and manipulations on the CAD file, in ASCII 
format, and translates it into a solver accept- 
able input format. The objects’ electromagnetic 
properties are included at this stage. 


3.4 The Post Processor Unit 

Under the supervision of a managing module in- 
terfacing the solver output database with the 
user control screen, the Post Processor Unit per- 
forms the following tasks: extraction of the time 
series of the received e.m. field, field weighting 
by the preselected mobile terminal antenna pat- 
tern, computation of the narrowband statistical 
functions and wideband channel parameters. It 
is important to stress that the user has the op- 
portunity to select his mobile receive antenna 
within a large set of predefined radiation pat- 
terns; a computed weight function can also be 
user defined to take into account the effect of 
the vehicle roof. In Fig. 4, the time series for 
the test case presented in Fig. 3 is reported. 

The set of statistical narrowband analyses 
available to the user includes Probability and 
Cumulative Distribution Functions (PDF and 
CDF), Average Fade Durations (AFD), Level 
Crossing Rates (LCR), Distribution of Fades 


and Connections, Time-share of Fades and Con- 
nections. 



T [SEC.] 

Fig . 4 Time scries (segment A-B of the 
trajectory) 


The list of wideband channel parameters is also 
quite exhaustive: Doppler and power delay pro- 
files, channel scattering function, delay spread 
and coherence bandwidth. Through the use of 
these routines the LMS system engineer can ef- 
fectively achieve a very detailed characterisation 
of the narrow and wideband channel, ([6]-[7]). 
Fig. 5 reports, on the base of the previous ur- 
ban layout, the corresponding PDF and CDF; 
it is clearly visible the multipath effect due to 
the strong reflecting buildings (perfect conduc- 
tivity) along the trajectory considered. 
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Fig. 5 PDF and CDF 
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In Fig. 6, the scattering profile of the urban 
channel computed in several observation points 
along the selected trajectory is presented. It is 
very interesting to observe the mutual interac- 
tion between different ray contributions gener- 
ated by the urban area elements. 



Fig. 6 The Scattering Channel Profile 

Finally, the cumulative distribution functions of 
the average delay and delay spread are given in 
Fig. 7. The segmented shape of the curves is due 
to the limited number of observation points con- 
sidered in our test case. It is fairly easy to notice 
that the maximum average delay is around 76 ns 
while the delay spread is always less than 40 ns. 

4 The hardware and soft- 
ware platforms 

The LMS prediction tool for urban environ- 
ments has been presently designed to run on 
PCs with MS-DOS operating system. The re- 
quirements for the hardware platform are an 
80486 processor, 8 Mbyte of memory and 200 
Mbyte of hard disk, a VGA or better EGA 
graphic card. FORTRAN has been used to de- 
velop all the model units but the Man Machine 
Interface running under Windows environment 


through Microsoft Visual Basic. The user can 
visualize and interact with the Urban Modeler 
via a DesignCAD-3D tool, incorporating a large 
set of commands and macros. 




Fig. 7 Average Delay and Delay Spread CDFs 

5 Conclusions 

A brief summary and description of the main 
features of a GTD-based LMS prediction tool 
for urban environments has been presented in 
this paper. The model can be effectively used to 
simulate and estimate the behaviour of the nar- 
row and wideband channel for any given mobile 
satellite system in built-up areas. The user can 
define his own urban layout through a Design- 
CAD tool inserting buildings, tunnels, phone 
booths, parked vehicles and many other typi- 
cal urban items. The model itself is suitable 
for simulating moving vehicles at given speed as 
well as personal communication network users 
equipped with hand-held terminals. The set of 
wideband parameters and statistical functions 
available to the user, being a propagation engi- 
neer or an LMS system planner, is fairly com- 
prehensive. 

In the next future, the LMS prediction tool 
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will be upgraded insofar as higher orders of elec- 
tromagnetic contributions and the transmitter 
located inside the urban layout will be taken 
into account and implemented. With proper 
but minor modifications of the Ray Tracer and 
the GTD solver and with the consequent up- 
date of the MMI, the model itself will be then 
able to estimate channel parameters also for in- 
building and cellular radio communication net- 
works. In addition to the afore mentioned im- 
provements, the LMS prediction tool will also 
undergo fairly extensive and comprehensive vali- 
dation campaigns in the ESA Compact Antenna 
Test Range, on a scaled model of urban environ- 
ment, and using actual experimental data. 
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